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StMtART 

An approximate method, for computing thermal stresses In a thin 
plate of variable thickness is presented.. The temperature and. the thick- 
ness can both vary chardwise and. s panwise, and the effect of the free 
end Is automatically included. The method makes use of polynomial ap- 
proximations for the stress function to reduce the partial differential 
equation to a set of ordinary differential equations. This results in 
satisfying the differential equation everywhere spanwise and at a finite 
number of stations chordwlse. The boundary conditions are everywhere 
satisfied. 

Several examples of the method are presented in detail, and curves 
showing the effects of chordwlse and spanwise temperature distribution 
and chordwlse and spanwise thickness variation are shown far several 
cases. It is indicated from these few examples that the effects of 
relatively large thickness variation and spanwise temperature variation 
are not of major importance for a plate with a free end. 


INTRODUCTION 

Jet engines, high-speed airplanes and missiles, and nuclear power- 
plants are examples of modern devices in which large temperature gra- 
dients exist. Such temperature gradients can produce large thermal 
stresses which, by themselves or in conjunction with stresses produced 
by various external loads, can cause serious component failures. Thus 
jet-engine turbine blades possess large temperature gradients during 
operation which cause internal forces that superimposed upon the centri- 
fugal and aerodynamic loads can contribute to blade failure. During 
acceleration, temperature gradients occur in wings of high-speed aircraft 
because of aerodynamic heating. Serious thermal stresses may thus be 
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produced which, must he of concern to the designer. ,These and other new 
applications have caused a renewed and urgent interest in the relatively 
old subject of thermal stresses. 

The geometries that exist in these various applications are usually 
very complicated and difficult to treat exactly. It is usual, there- 
fore, to approximate the various configurations by simple bodies such 
as plates, cylinders, and spheres. Thus, the turbine blade and the 
supersonic wing are approximated in this investigation by the t hin plate 
of variable thickness. 

The thermal stresses in a thin plate are treated in detail in ref- 
erence 1. However, the plate is assumed to be of constant thickness, 
the temperature is assumed to vary chordwlse only, and the end effects 
are neglected, use being made of Saint Tenant's principle. All these 
assumptions are invalid in many practical cases. In a turbine blade, 
for example, there is a large variation both in temperature and in 
thickness in both chordwlse and spaavise directions. In addition, the 
length-to-chord ratio is often less than 2, and the end effects may 
therefore be Important over a large part of the blade. It is suggested 
in reference 1 (p. 400) that the end effects can be calculated using 
the strain energy method. 

In order to avoid some of these assumptions, an approximate solution 
to the equations of elasticity is given in reference 2. A solution is 
obtained far the bl harmonic equation far the stress function in terms 
of an Infinite series. However, only two terms of the infinite series 
are used, and the boundary conditions are not exactly satisfied but only 
in an average fashion over the surface. The results at the free end of 
the plate are compared with those obtained by the energy method, and 
it is indicated that the method of reference 2, although mare labarous, 
has inherently higher accuracy. The calculations made in reference 2 
are for a uniform- thicknes s plate with chordwlse temperature variation 
only. It is indicated that, although this method, as well as the energy 
method, could be adapted to a variable-thickness plate, the amount of 
labor involved would become prohibitive and should not be attempted 
without the aid of a high-speed computing machine. 

The equations for calculating thermal stresses in a supersonic wing 
with chordwlse variations in thickness and temperature are presented in 
reference 3. The method of reference 1 (p. 399) is used in deriving 
these equations, and again the results as stated therein are not valid 
at a distance less than a chord length from the tip. 
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The object of the present paper is to present an approximate method 
for calculating accurately the thermal stresses In a thin plate with 
both chordwlse and spanwlse variations In temperature and chardvlse but 
restricted spanwlse variations in thickness. Numerical examples of the 
effects of the free end and of thickness variations will be presented. 

The method used Is essentially a collocation procedure applied to a 
partial differential equation and Is used successfully in reference 4 
for calculating transient thermal stresses. In this procedure, the 
differential equation is satisfied everywhere in x but at only a finite 
number of values of y. Use is made of a polynomial approximation for 
the stress function to reduce the biharmonic equation to a set of or- 
dinary linear differential equations which can be readily solved. The 
extension to other problems involving the bl harmonic or Laplace equation 
can be made. 


ANALYSIS 
Stress Function 

Consider a thin flat plate of variable thickness, as shown in fig- 
ure 1. The origin of coordinates is taken at the middle of the free 
edge, and the midplane of the plate is assumed to He in the xy-plane. 
The width of the plate runs from y =» -1 to y «■ 1; and the length, in 
the x -direction, can be either finite or infinite. The variable thick- 
ness is a function of both x and y and is assumed to have the fol- 
lowing form: 


h = h^xjhgfy) 

The function ^(y) is any continuous function of y having a continu- 
ous first derivative, and the function h-^x) is a function of x hav- 
ing the form h-^(x) «= h^e 1 ™. (Symbols are defined in appendix A) . For 

m equal to zero, the thickness is constant in the x-dlrectlon. At any 
instant of time, the plate is assumed to have a temperature distribution 
varying with x and y but independent of z. Thus 

T * T(x,y) 

It is assumed that the plate is sufficiently thin so that a state 
of plane stress exists, the only significant stresses being a x , a y , 

and n^., which are independent of z. The assumptions that 

a z “ T yz = \z ° 0 and that cr x , cr y , and are independent of z 

lead to an inconsistency in that all the compatibility conditions are 
not necessarily satisfied. It can be shown (see, e.g., ref. 1, p. 24l), 
however, that the error is proportional to z^ and is therefore small 
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for thin plates. It Is further assumed that any curvature of the mid- 
plane of the plate can he neglected. 


Under these conditions it is shown in appendix B that the elastic 
stresses a x , o^., and x^ are derivable from, a stress function <P 

satisfying the following differential equation: 


^1***“)* (i + v ) (e) 


4> + 




where the subscripts indicate differentiation with respect to that vari- 
able, and v 2 is the Laplnclan operator 


V 2 


= a 2 

i? dy^ 


a- « r- 


The stresses are given In terms of the 

1 

_ la 2 ? 


stress function as follows 
*\ 

► 


T 1 ^ Z f 

h chTSy J 


(2) 

or 

(B5) 


The boundary conditions are ones of no normal or shear forces on 
the surface of the plate. These conditions lead to the following equa- 
tions that must be satisfied at the edges of the plate: 


at x «= o 


at y =» ±1 


h0_, 


hT. 


xy 



- d 2 ? 

dx dy 



hoi 


hx. 


xy 


d 2 ? 

a?" 

. d 2 ? 

dx dy 



( 3 ) 
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For an infinite plate the stresses at Infinity must be finite. For a 
finite plate with both ends free, the same boundary conditions exist at 
both ends. 

For a fixed edge the boundary conditions are given In terms of 
displacements. This necessitates setting up the differential equations 
In terms of displacements rather than stresses. 


Polynomial Approximations 


Equation (l) will be solved approximately by a collocation proce- 
dure whereby the differential equation Is satisfied everywhere in x 
but at only a finite number of values of y. Thus, as Indicated in 
figure 1, n stations are taken along y. The y-coordlnates of these 
stations are y y 2 , . ... y Q . The stress function q> is then assumed 

to have the following farm: 




(s) 


where Pj (y) Is a polynomial in y associated with the J ttl station 
and satisfying the following conditions; 


V.VUVLL UAvUB • 

" 1 \ 
r t ) - o AjiJ 


(6) 


Since will. In general, be a nonlinear function of x, the fol- 
lowing conditions must be satisfied in order to ensure that equations 
(4) hold far any <Pjj 


Pjtti). 


:} 


(7) 


Polynomials having these properties can readily be obtained. Thus 

yj jf±i 


& i7 ' ' y i> 


( 8 ) 
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where n Is the product for a.11 values of 1 except in], Equa- 
led 

tlou (8) satisfies equations (6) since Pj » 0 at y ■» y ^ and. Pj «• 1 

at y o yj. It also satisfies equations (7) since Pj and Pj are 

hoth zero at y = ±1. Equations (8) and (5) are new substituted into 
equation (l) , and the equation is then evaluated at each of the n 
stations. This results In a set of simultaneous fourth-order ordinary 
differential equations in (x) of the following form: 

C a ij*J n + *1^3" + Cij’S + *± 3 *$ + “ -h(x J y 1 )V z |w(x,y 1 )| 


1 = 1, 2, . . . , n 


(9) 


where 

- p 3&i> 

h 5x 


“ ' sr h i p j( y i ) 

/ I* ^1 ^ . 

, o <2EP3 - ^ + pJ4 (Of - i S& - m* + y. ^li 

iJ \ J 11 h ^x 2 h 2 ^/ h dy^lj 3 

- {»; - K3T + v ® 

a lJ“ 2 [-#S + (l + v)p j(^^S'E^)] J-yi 

" 2 &S +(i+v,p iS| wl 

■ - {t - - T | ■ * > n [£ (If - k$ - 5 (i 2 ♦ e gj} 

l p r - *r 3 + ^BsT * 3 - *0 * v S} W1 


>(io) 
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It Is to be noted that because of equations (6) an terms in equa- 
tions (10) which are multiplied by Pj are zero except when J = i, 

in which case Pj (yj ) = 1. 


Special Cases of Thickness Variation 


The coefficients defined in equations (10) are, in general, func- 
tions of x if the thickness varies s panwise. The solution of equa- 
tion (9) can became very laborious under these conditions. If, however, 
h^(x) is chosen to be an exponential li^e 312 , then the coefficients are 

all constant. Equation (9) can then be solved subject to the boundary 
conditions (eqs. (3)) by any of the standard methods for solving ordi- 
nary linear differential equations with constant coefficients. The so- 
lution will be given in terms of exponentials. By the proper choice of 
m, h^(x) can be made to give a reasonable approximation to the variation 

of thickness along the span of a turbine blade, particularly in the 
vicinity of the free edge. If there is no variation in thickness, span- 
wise, then m is taken as zero. 


If the thickness varies only in the chordwlse direction, b. . and 


ft 


and all terms in 
then becomes 


ij 

Cjj and e^j containing h^ vanish. Equation 


(“13*3” + c i 3 ,, 3 + *13*3) ■ -fcCy^lwfoyi)] 


(11) 

i*=lj 2 j • » • 


If the thickness is constant throughout the plate, equation (l) 
becomes 


V 4 <p « -V 2 (EaT) 


( 12 ) 


where 



a 4 , - a 4 1 a 4 

a* 4 h? a^ 8y 4 


is the blharmonlc operator, and h is deleted from equations ( 2 ), which 
define the stress function q>. Equation (9) reduces to 



« 


(7^ 


+ 2Pj(y.i>3 + 


(13) 
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Special Cases of Temperature Distribution 

If the thickness and the temperature are even functions in y (i.e., 
symmetric about y «* 0), then only that half of the plate between 
y« 0 and y ** 1 need be considered. This will reduce the number of 
stations required for a given accuracy. The polynomials P . can be 
taken as even polynomials; thus J 

p jW - - & * 3 ^ tu) 

J 

Similar ly, if the thickness is an even function of y and the tempera- 
ture is an odd function of y (antisymmetric ) , the polynomials are taken 
as odd polynomials; thus 

p jM - % ~ ^2 i5j TCy2 ‘ V* 7j(5 J ' ^ y J * (1S) 

Again only that half of the plate between y » 0 and y » 1 need be 
considered. Since any tenqserature distribution can be split up into 
the sum of an even and an odd function, it is possible to solve the 
problem of an arbitrary temperature distribution in two stepB, namely, 
solving for the stresses far an even distribution and far an odd dis- 
tribution, and adding the two. This will generally result in some re- 
duction in labor, since the work Involved in solving two problems each 
with n/2 stations is usually less then solving one problem with n 
stations. It is to be noted, however, that this procedure cannot be 
followed if h la not an even function of y. In that case the general 
polynomial (eq. (8)) must be used and stations taken ranging from y ■» 1 
to y - -1. 


EXAMPLES 

Uniform Thickness, Parabolic Chardwlse Temperature Distribution 

As a first example, consider a send -Infinite thin plate of constant 
thickness with a temperature distribution independent of x and given 
by 

T = (y 2 - |)T q 

Since this temperature distribution has zero mean and zero moment about 
the x-axis, the stress at a distance far from the free end is given by 


4098 
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the method of reference 1 (p. 401 ) as 


cr = -Eaff 

Jfc 


^3 - fyo 


a = t s 0 

y zy 

The only problem, therefore. Is to obtain the stresses near the free 
end. Since the thickness is constant and the temperature distribution 
is symmetric, equations (13) and (14) vill be used. Two stations are 
chosen at y^ = l/4 and y 2 = 3/4. The polynomials then became 

(y 2 - 1 ) 2 (y 2 - jg) 512/2 . 2/2 9_\ 

Pl ° 72 ^WT JM " 225^ " ^ V “ 16/ 

\16 " 7 \16 “ 16/ 

(y 2 - 1 ) 2 fy 2 - 


w V 16J 

/JL . A*(± . A 
\16 7 \16 16 J 


§*&* - ^ - £) 


Evaluating P-, and P g and their derivatives at y, =■ 1/4 and 
at y g =■ 5/4 gives c 


P-lU/ 4) = 
P£(l/4) ® 

p™ (1/4) 

P 2 (l/4) . 
P"(!/4) * 


1 

-5.564 
» 88.75 
0 

8.571 


^(3/4) «B 0 
P£( 3 / 4 ) B 8.089 
P™ (5/4) « -320.8 
P 2 (3/4) - 1 
Po ( 3 / 4 ) » -22.78 


Pg M (l/4) « -282.1 Pg" (3/4) = 1599 
Substituting now into equation (15) gives 

<P™ - 11.13<p£ + 88.75^ + 17.14<pg - 282.1^ » -2EaT 0 

16 . 18<p£ - 320. 0P 1 + <P 2 n - 45.5502 + ^"^Z “ -2BxT 0 
The particular solutions are readily obtained as 


«P lp b -0.07324EoTq 
<P 2p ■ -0.01595 EocTq 
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The homogeneous solution Is obtained "by the usual exponential sub- 
stitution. The complete solution Is then 


<Pl = 



0.07324EoTq 





0 . 01595Ea3? 0 


■where the terms X^ are the roots of the determinants! equation and are 
equal to 


^ « -2.120 + 1.1171 

*2 * *1 
*5 * 

*6 * *5 


X 3 » -5.682 + 2.6811 

X 4 ” *3 
X 7 » -Xg 


where the bar signifies the complex conjugate. 


The term Is given In terms of as follows: 

xj - 11.13x5 + 88.75 
Bjj.- - -t _J5 

17.14XJ - 282.1 


In order for the stresses to remain finite as x approaches infinity, 
Ag, Ag, Ay, and Ag must vanish. This leaves for and q> 2 


«p 1= . 



0.07324BoT 0 


xj - 11.13X5 + 88.75 Xjx 

*2 - •“ > — 2 ~— A fc e - 0.01595 EccTq 

fel 17 . 14X^ - 282.1 


In order to determine the values of the A k 's, use is now made of 

the boundary conditions of equations (3) . It is to be noted that all 
other boundary conditions are already satisfied in the choice of the P 
function. In order to satisfy equations (3), it Is necessary that 

<Pl(0) = <P£(0) - q> 2 (0) a ¥£(0) - 0 


/ 
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co 



since Six Is not necessarily zero at x = 0. Thus, four equations 

by 2 

result for the four A^’s. The solution of these four equations gives 

A x = (0.03625 - 0. 071141 )BccT 0 Ag = A^ 

Aj » (0.0003671 + 0. 00019261 )MP q A 4 - A^ 

After adding all the complex conjugates, the Imaginary parts drop 
out and the final solution is 

- 2e" 2,12Qjc [0. 03625 cos(l.ll7x) + 0.07H4 sin(l.U7x) JBoTq + 

2 e "5.682X[ 0>0003671 cos(2.681x) - 0.0001926 sin(2.681x) ]EcxT 0 - 
0.07324EocTq 

(Pg « 2e" 2 * 120x [0.008941 cos(l.ll7x) + 0.01249 sin(l.H7x)]SiT 0 - 

2e -5.682xj- o ooog662 COB (2.681 x) + 0.0001790 sin(2.681x) ]EoT 0 - 
0.01595EcxT- 

The stress function Is given hy 

<p» P-l^ + P 2 <P 2 

and the stresses are 
ff x ° P 1*L + 

«T “ p l*£ + p 2 *1 

- - S&* - + A+ - rh 



-P ’«P 1 - P 'Vg 

§§(^ -V - $)* i '?te 
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It Is to be noted that at x ** « 


Therefore, 


(pj_ «* -0.07324 Bo3?q 
C fVj a -0.01595Eaa?Q 

«Pi - = <p^ *= <Pg = 0 


*7 

(-y 2 + 0.333)Bcdr 0 


which are the exact stresses at x a «• . The particular solution of the 
differential equation gives the stresses far from the end. The homo- 
geneous solution gives the correction due to the free end and dies out 
rapidly away from this end. 


The preceding stresses are plotted in figures 2 and 3. These results 
will he discussed in BESQHES AND DISCUSSION. 


The values of the X's, given previously, can he used far any 
constant-thickness two-station symmetric problem, since the X's are 
independent of the temperature distribution. Far any other temperature 
distribution new values far the particular solutions to the differential 
equation must be obtained, the homogeneous solution remaining the same . 
New values of the constants and can then be determined from 
the boundary conditions. 

Far a three-station solution of the same problem ^stations at 
7 ± 01 -g, yg *= and y^ a -gj, the following values can be used: 

PiCy) * 7.14l(y2 - l^y 2 - i^y 2 - ||) 

P 2 <7) - -1B.OO& - !) 2 (/ - iV/ - §§) 

P 5 ( y ) - 36.15 (y 2 - U^y 2 - ^(y 2 - i) 

\ a -2.10568 - 1.125091 Xg » X x 

Xg a -5.42265 - 1.278261 X 4 a \g 

Xg a -9.23538 - 3.371451 Xg - Xg 
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The equations for the <Pj functions In this case are 


ip 


A . v 


1 - 2^ V + ''ip 



* 2 - ? Vk e + <p 2p 


V V 

** “ 2./ Vt e + % 


where the tPj^ 'b are the particular solutions for the given temperature 
distributions and the K^ and values are as follows: 


% => 0.5637 - 0.039471 % 
Kg a -1.236 - 0.25411 K 4 = % 
Kg a -2.968 - 1.0611 Kg “ % 
K-l » 0.08273 - 0.018021 H g » 5^ 
K 3 « -0.4420 + 0.040441 = 5g 
% « 2.019 + 4.4971 Hg a S g 


the A^'s can be obtained from the boundary conditions; that la, 

«Pj a <pj o 0 at x a 0. 

Variable Chordwlse Thickness, Parabolic 

Chordylse Temperature Distribution 

As a second example, the plate Is assumed to have a variable thick- 
ness In the chordwlse direction given by 

h a hg - (1 - 0.9y^)hQ 
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*T 1 hf> temperature distribution Is taken the same as before: 

1 - if - iK 

Again two stations will be taken at y-j_ » l/4 and at y 2 “ 5/4. Since 

conditions of symmetry exist, the same polynomials as used In the first 
example will be used here. H'hw derivatives of these polyp.om1.alB eval- 
uated at the two stations have previously been listed. Since the thick- 
ness Is no longer constant, equation (15) cannot be used. Instead, 
equation (ll) is used In this case. The coefficients a^j, c^j, and 

are evaluated using equations (10) . Equation (ll) then becomes 

cp™ - 13.81<p£ + 104. ^ + 21.52q>2 - 366.5q>g » -l.SSShoBtxTo 

12.61<p£ - 226.4<P 1 + Vg' - 72.22<pg + 1947q> 2 - -0.9875h 0 BaT 0 

These equations can be solved In the same way as for the uniform- 
thickness plate. The final results far <p^ and <Pg are 

2e -2.569xj- 0>01644 coa (i,;i 59 x) + 0.03761 sin(l.l59x)] + 

2 e -6.386X[ 0>()002130 C06 ( 1#741x ) + 0.000002496 sin(l.741x)] - 
0.05330 

2e -2.569x[ 0 QQgggg cOS (x.l59x) + 0.003692 sln(l.l59x)] - 

2 e “ 6 * 38ex [ 0.0004389 cos(l.741x) + 0.0001872 sln(l.741x)] - 
0.004381 


h 6 Ba!C 0 


9 2 


IIqBxTq 


The stresses are then given by 
ff x " h P 1 9 1 + h P 2 9 2 


512 

225 


( 30 , 4 . 2 - 2 .-) + 

(l-O.Sy 2 ) ^ 49 (l-0.9y^ ' 
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“y - £ F l*i * j 

= _ 512 (y 2 - - jsL S12 fa 2 - X ) 2 (^ - mL 

~~ a - o vi ^ + « (i . oV) ^ 


225 


V - - E p i*i - 5 P 2*2 


512 

225 


(earS-^uft) 512 (gy 5 - fy 2 


(1 - 0.9y 2 ) ^ 49 


(1 - 0.9y ) 


At x o 


a v s* *c = 0 

y xy 


(5Qy 4 - 30.75y2 + 4.25) _ 


~~ - 0.07578 — «- 

^0 (1 - 0.9^) 


0 . 0*578 Ssedi^jsd **■*). 

(1 - o.gy 2 ) 


- 0.219 - y 2 

The exact solution far this case at x « » Is given by (ref. 3) 

A * hay 
T . J- 1 T 0 


x 

BaT, 


o t ° r hay 


Substituting the expressions far T and h gives 


EaT, 


-y2 + i + 


0.219 - y 2 


i jC (»* ‘ s) (1 - °- 9,2)a3r 


X 0 -"- 


Thus the two-station solution again gives the same answer as the 
exact solution at x equal Infinity. 

The stresses x xy’ 8114 a x for this variable-thickness case are 

plotted In figure 4. These results will be discussed In HBSUIT8 AND 
DISCUSSION. 
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Uniform Thickness, Honsymmetrlc Chordwise 
Temperature Distribution 

As a final example, consider the case of a uni form- thickness plate 
with a nonsymmetric chordwise temperature distribution given by 

i - (j 3 + y 2 - s)t 0 

Since this temperature distribution is neither even nor odd., the general 
polynomial of equation (8) will be used with three stations taken at 
y-^ a -2/3, =* 0, and y^ » 2/3. The polynomials then are 

- SS 4 - D&* - 1)2 


K* + D( y ■!)&*- 1)2 


r(y + | - 1) 2 


When these polynomials and their required derivatives are evaluated at 
the three stations and the results substituted into equation (13), the 
following three equations are obtained: 

•p™ - 30.0ep£ + 797.0 <p l - 228. (XPg - 9.90Gp£ + 213.8<P 3 = 2BaT 0 

14.58(p£ - 175. 0^ + <p£* - 17.00q>g + 152. 0q>g + 14.58q>3 - 175.0q> 5 » -2BdT 0 

-9.90(XP£ + 408.2^ + lS.CXXPg - 228. 0q^ + <pg' - 30.06«p£ + 602.6<P 3 » -SEoTq 

These equations are solved giving <P^, «Pg, and <P 3 as functions of x. 

The stresses are now computed as before from equations (2) with h 
deleted. 

The previous problem could also have been solved by separating 
the temperature distribution into even and odd functions. Thus the 

stresses far the distribution T ■ ^y 2 - — have already been com- 
puted in the first example. It remains only to obtain the stresses due 
to the distribution T = j^Tq. This can be done accurately by using 

only two stations at y^ = 1/4 and yg *» 3/4 with the antisymmetric 
polynomials of equation (15) . The two stress distributions can then be 
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added to give the complete stress distribution. This has been done for 
this problem and the results are compared vlth the previous solution in 
figure 5. Figure 5 also Includes a four-station solution of this same 
problem using the nonsymmetrlcal temperature distribution. 


RESULTS AHD DISCUSSION 

In order to determine the relative accuracy of the method presented 
herein, a comparison is made in figure 2 between this method using one, 
two, and three stations and the energy and elasticity methods presented 
in reference 2. A uniform-thickness plate with the parabolic tempera- 
ture distribution of the first example Is used for comparison. The 
spanwise stress c x , the shear stress at the tip (x = 0), and the 

chordwlse stress cr at l/4 chard from the tip (x « 1 / 2 ) are plotted 

¥ 

far the three methods. It Is seen that far an practical purposes the 
method presented herein has converged using Just two stations, giving a 
far superior answer than the energy method and a slightly better answer 
than the elasticity method of reference 2. Even the one-station solution 
gives an approximate answer almost as good as the energy method and may 
be useful far many practical engineering problems. 

The effect of the free end for this problem is briefly shown in 
figure 3(a), where a x and at the midchord (y * 0 ) are plotted 

against x. The chordwlse stress o y is large at the tip but has 

dropped to practically zero after 1 chord length from the tip. The 
spanwise stress o x starts at zero at the free end and rapidly rises 
to a constant value after 1 chord length from the tip. The shearing 
stress, although not shown here, follows the same pattern as cr y . It 

Is to be noted that an approximate simple check can be obtained far cl. 
from the fact that the integral under the curve must vanish. It is 
further seen from this curve that there is practically no difference 
between the two- and three-station solutions far this problem. The com- 
plete stress distributions far this case at various distances from the 
free end are plotted in figures 5(b) to (d). 

The effect of variation in chordwlse thickness Is shown In figure 
4. Again the temperature distribution is parabolic and the thickness 
variation Is given by h = ]1 q( 1 - 0.9y2), where hQ is the thickness 

at the m id ch ord. The edges therefore are one-tenth as thick as the 
middle. A comparison Is made between a two-station and a three-station 
calculation, and it Is seen that the two-station calculation gives good 
results except for the spanwise stress a x near the thin edge, where 

there is some difference between the two- and three-station solutions. 
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A comparison 'between the uniform-thickness plate and the variable- 
thickness plate in figure 4(d) indicates, as would he expected, that the 
effect of thinning out the edges is to raise the edge stress and lower 
the stress at the midchord. However, in spite of the large variation 
in thickness, the edge stress is raised by only about 15 percent far 
this case. 

The case of a nonsymmetric chordwise temperature distribution 
T =» + y^ - ^ (the third example) is shown in figure 5. Three 

methods were used for this example: a three-station solution; a four- 

station solution; and a solution obtained by splitting the temperature 
Into even and odd functions, using a two-station solution for each tem- 
perature distribution, and adding the results. It is seen from the 
figure that the sum of the two-station solutions gives the same results 
as that for the four-station solution. It is generally less work to 
perform two two-station calculations than one four-station calculation. 
Therefore, if the chordwise temperature distribution is not symmetric, 
it is best to divide it into even and odd parts and perform two separate 
calculations. This should be done if extreme accuracy is desired. 
However, it is seen from figure 5 that even three stations give suffi- 
cient accuracy for most engineering purposes; therefore, using more 
stations or dividing the problem Into two parts would generally seem 
unnecessary. ■ 

The effect of variations In both chordwise and spanwlBe temperatures 
Is shown in figure 6. The thickness is assumed constant and the tem- 
perature distribution Is given by T =■ T Q ^y^ - ^y(l + 0.3e _x ) . This 

would correspond to the temperature decreasing exponentially by about 
30 percent from the tip to the base of a turbine blade. This relatively 
large variation in spanwlse temperature affects only the chordwise 
stress o T , leaving the spanwlse stress cr x unaltered. Since the max- 
imum stress is generally the spanwlse stress a x , the effect of the 

spanwlse variation in temperature doesn't seem to be of great importance 
and can be neglected as a first approximation. 

As a final calculation, the effect of spanwlse variation of thick- 
ness was briefly Investigated. The temperature distribution was taken 
to be parabolic and the thickness was assumed to be given by 

h a hi = hge® ,3x . This is a reasonable type of spanwlse thickness vari- 
ation for turbine blades. The results are shown in figure 7. A com- 
parison of figure 7 with figure 3 shows that this variation in spanwlse 
thickness produces no large changes in the stress distribution. 
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Although the effects of thickness variation are shown to be not very 
great In these examples , this conclusion may not be valid for extreme 
thickness variation. 

® It Is to be noted that, although the cases considered herein are 

<3 all for a plate with a free end, any other edge conditions can be 

treated in the same way, provided the boundary conditions are adequately 
specified. 


ccncLusiars 

An approximate method was presented for calculating the elastic 
thermal stresses In a thin plate of variable thickness. From the 
examples shown. It would seem that a two- or three-station solution 
should suffice far most practical engineering problems. 

Calculations on the effects of spanwlse variations in temperature 
and spanwlse and chordwlse variations in thickness Indicated that for 
the examples chosen these variations did not have any large effects on 
the stresses. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 16, 1956 
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APPHHXEX A 
SXMBQL8 


■u ^ 
b u 


C iJ > 


■u 


h 


constants cf integration 


coefficients in differential equation for stress function 
(eq. (9)) 


Young's modulus of elasticity 
thickness of plate, h « ^(xJhgCy) 


*0 

\(x) 

^(y) 

*k 


X 


y 


reference thickness 

function expressing sparwise thickness variation 
function expressing chardsise thickness variation 
constant 

coefficient in equation h^x) *» 
constant 

polynomial associated vith station, function of y 
temperature, function of x and y 
reference temperature 

spanwise coordinate measured from free edge, h al f chard lengths 

chardwise coor dinat e measured from centroid, y-axis coincident 
vith axis of minimum moment of Inertia 


z 


thickness coor dinat e measured from center, z-axis coincident 
vith axis of maxi mum moment of inertia 
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a coefficient of linear thermal expansion 

T™ shearing strain 

Ajr 

s x strain In x -direction 

e strain In y-dlrection 

J 

v Poisson ratio 

II product for all values of 1 except 1 » J 

1 4 J 

<r x normal stress In x -direction 

cr normal stress in y-dlrection 

shear stress 

<p stress function, function of x and y 

q>j stress function associated with J til station, function of x 

<Pjp particular solution to differential equation 

Laplacian operator, 2 — + 2L_ 

3x 2 By 2 

V 4 hlharmonlc operator , + 2 — 

ax 4 ax 2 ay 2 ay 4 


Subscripts : 

1, J summation or multiplication dummy indices or refer to the 

1 th or J t]l station 

x,y partial differentiation with respect to that subscript, except 

where otherwise defined 

Superscripts : 

1 Indicates ordinary derivatives 

” indicates complex conjugates 
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APPMDIX B 


DIFFERENTIAL EQUATIC H FOR STRESS FUNCTION 
The stress -strain relations for the plane stress problem are 

E8 Sc “ “ vo y + 331 

a„ - va„ + EotT 


'xy 


~7 x 
2(1 + v) 


(Bl) 


E 


''xy 


The compatibility equation Is 

a!sx + ?!fz 
a * 2 a * 2 

Substituting equations (Bl) Into equation (B2) gives 

V2(a x + a y ) . -V^(BxT) + (l + v)Q-jj£ 


(B2) 



(B3) 


The equations for the equilibrium of forces in the plane of the plate 
are (ref. 5, p. 297) 


d(h<x x ) d(lrO 

~ 3F“ + ~=P £ ~ 


W 


d(bQy) d(hT_) 

If a stress function 9 Is defined as follows: 

‘ 1 B 2 <P 
or o — ■ - ■ ■ 

o . 1 8 2 * 

* *5? 

nr - - 1 ^ 2< P 


(B4) 


(B5) 


then equations (B4) are automatically satisfied. Substituting equations 
(B5) into equation (B3) gives the equation to be solved for the stress 
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function: 


V 2 <p) - -V 2 (BoI) + (1 + 4 <Pyy)xx + (e <boc)yy - 2(e <ky)xy| 

m 

»2| V 2 *] - Cl + ^ + (sL'te - 2 (e) «»0 - -^CW) (B6) 
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(a) Comparison of stresses at midchord for two- and three-station solutions. 


Figure 3. - Stresses in plate of uniform thickness far parabolic temperature 
distribution. T - T^y 2 - . 


* 


c 


l 



Dimensionless chordwlse stress 
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.0 

Chardwlse distance measured from centerline, y, 
semi chorda 


(c) Spaxnrlse stress; two-station solution. 

Figure 3. - Continued. Stresses In plate of unif orm 
thickness for parabolic temperature distribution. 
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(d) Shear stress; two-station solution. 

Figure 5. - Concluded. Stresses In plate of uniform 
thickness far parabolic temperature distribution. 

T - T^y 2 - |). 
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J f . , l. , A. — — | . | . ■■ —I— — M -l— — ■■ ■ I I I 

* 0 .2 .4 .6 .8 1.0 

ChordwlBe distance measured from centerline, y, 

semi chor da 


(a) Chardsrt.se stress. 

Figure 4. - Stresses In plate of variable chordvrt.se 

thickness, h » 1]q(1 - O.Sy 2 ); T *» TgCy 2 ~ ?)• 

o 


4098 





.6 .8 

Chardvise distance measured frc* centerline, y, 
semi chords 

(b) Spamrise stress. 


Figure 4. - Continued. Stresses In plate of verl 
able chardwise thickness, h - hgfl - 0.9y 2 ); 

T-Tofr 2 -! 5 * 
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Stations 
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» 
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— *2 


Chardwlse distance measured, from centerline, y, 
semi chords 


(c) Shear stress. 


Figure 4. - Continued. Stresses in plate of vari- 
able chordwise thickness, h = 1 iq( 1 - 0.9y^)j 

t « TqCj 2 - |) . 
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ChardwlBe distance measured from centerline, y, 
semi chords 


(t>) Chordwise stress. 

Figure 6. - Continued. Stresses in plate of uniform 
thickness for parabolic temperature distribution. 

T - ToCy 2 - i) (1 + 0.3e" x ) . 
o 


OB 
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Dimensionless spanwlse stress 
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semi chords 


Ob) Spanwlse stress. 


Figure 7. - Continued. Stresses In plate of vari- 
able thickness . h a T a - i) . 


Dimensionless shear 



Chardwise distance measured from centerline, y, 

semi chords 


(c) Shear stress. 

Figure 7. - Concluded. Stresses In plate of vari 
able thiclmess. h a hoe° ,3x ; T = Tgfy 2 - 
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